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1 Introduction 



The study of rescattering effects in heavy flavour production in hadron-nucleus and 
nucleus-nucleus collisions is of particular interest for the interpretation of J/ip suppres- 
sion found experimentally at CERN It is well known that an important part of this 
suppression is due to the rescattering of the pre-resonant cc system with the nucleons of 
the colliding nuclei. This phenomenon is known as nuclear absorption. To describe it 
the following probabilistic formula [@, H is currently applied 



n'^ — n'^ 



/r+oo / r+oo 

dh dz pA{b, z) exp -cr^^^ A / dz' pA{b, z'\ 
J —oo \ Jz 

J d% [l - exp (-at^ATAib) 



— (JpNy^abs 

In this formula pA and Ta are the total and transverse nuclear densities respectively, 
^pA{N) inclusive J/^j cross-section on A{N) and a^^^ is the absorptive {cc)-N 

cross-section (i.e. corresponding to final states without J/ip). For open heavy flavour 
production aabs = and one obtains an behaviour in agreement with experiment . 
For charmonium and bottonium Cabs 7^ and one obtains a behaviour A°' with a < 1 
- also in agreement with experiment f!], |^. Eq. (1) has a probabilistic interpretation 
with a clear longitudinal ordering in z: in the first interaction at z the heavy system 
is produced and in successive ones a.t z' > z it rescatters with nucleons along its path. 
However, there is a fundamental change in the physical picture as one goes from low 
energies, for which the probabilistic formula (1) is derived, to high energies. Instead of 
successive interactions of the projectile with nucleons of the target, one has simultaneous 
interactions of particles into which the projectile has split. Thus the nuclear interaction 
responsible for the heavy particle production may occur after the interactions responsible 
for its absorption in nuclear matter. To illustrate this somewhat paradoxical point, 
compare Figs, la and lb, in which the heavy particle production and its subsequent 



interactions is presented as seen in the lab. system. The time goes from left to right. 
It is clear that in both cases the heavy particle is in fact created long before the actual 
interactions with the nucleus. In Fig. la the times of these interactions correspond to 
the common idea of production first, absorption after. However in Fig. lb the order 
of interactions is the inverse one: the interaction of the (already formed) heavy particle 
precedes the one which is responsible for its formation. 

In what follows we are going to study the changes in Eq. (1) resulting from this 
change in the physical picture at high energy. We shall derive an expression valid at all 
energies which exactly coincides with (1) in the low energy limit. 

Our expression contains not only the rescattering of the heavy system but also the 
rescattering of light particles (gluons and light quarks) . The latter can be interpreted as 
shadowing corrections to the nucleus structure function and cancel out in the low energy 
limit. Also our formula splits into pieces with no interactions of light particles, which 
we call the internal components and which are shown to be small in the central rapidity 
region, and a piece with at least one light particle interaction, which we call external and 
which corresponds to the usual production mechanisms such as gluon-gluon fusion. 

At asymptotic energies the rescattering of the heavy flavour in the external component 
has a form different from (1). The change consists in the substitution 

ATA{b)exp(-hATA{b)y (2) 

where a is the cc — N total cross-section. It is most interesting that for a ~ (Jabs, which is 
the case for charmonium or bottonium production, the first rescattering correction (i.e. 
the term proportional to a) is the same in the two cases. Since a is small the change with 
increasing energy in the absorptive corrections due to the rescattering of the heavy system 
will be comparatively small. As a consequence, the probabilistic expression remains 
approximately valid even at LHC energies (if the shadowing of the nuclear structure 
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-cr: 



ATAib 



functions is neglected). 

The situation is completely different for open charm or bottom production. In this 
case aabs = and, as we said before, (1) leads to A^. However with increasing energy there 
are shadowing corrections resulting from the rescattering of the heavy system (despite 
o'abs = 0). At high energies {y/s > 200GeV) they are practically identical to those for 
the J/ip- This is a main prediction of our approach. 

On top of these effects due to the rescattering of the heavy system, there are, of 
course, corrections due to the shadowing in the nucleus structure functions. The latter 
are present at high enough energy for all systems including Drell-Yan pair production, 
but vanish in the low energy limit. 

The paper is organized as follows. In Section 2 we derive the expressions valid at 
asymptotic energies which show the change (2) in the absorptive series for the external 
component. These formulae are in fact applicable only at extremely high energies. In- 
deed, due to the presence of the heavy system there are finite energy corrections which 
are important up to energies of the order M^i?^/x+, where M is the mass of the heavy 
system, x+ its light cone momentum fraction and Ra the nuclear radius |^. These 
finite energy corrections are explicitly computed in Section 3. In this Section we also 
show that, in the low energy limit, our formulae coincide exactly with the probabilistic 
formula (1). Section 4 contains our numerical results. Our conclusions are summarized 
in Section 5. 

2 Heavy particle production at asymptotic energies 

In this Section we study the inclusive cross-sections for the production of a heavy flavour 
system in a high-energy collision, including the rescattering of the heavy system. We 
start with pA collisions. The total cross-section for production of the heavy system 
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can then be represented by the diagram shown in Fig. 2 p[. It has a clear meaning 
in the laboratory reference system: at some point the heavy system is born from light 
partons (gluons), whereafter both heavy and light particles scatter on the nucleus. In the 
probabilistic treatment leading to Eq. (1) there is a time ordering: first an interaction 
of a light particle occurs, in which the heavy system is created, and only afterwards the 
latter interacts with the nucleus. As explained in the Introduction this time ordering 
does not exist at high energies where all interactions of light and heavy particles are 
simultaneous. We shall treat interactions with the nucleus of both the light particle and 
heavy system in the eikonal approximation. We also take all participating particles as 
scalar, for simphcity. 

2.1 The external contribution 

Using an eikonal model for the multiple scattering (both of light and heavy systems) in 
Fig. 2b, we obtain for the external contribution to the inclusive cross-section of a heavy 
system, integrated over the transverse momenta: 

where 

<^A = -[(^A{a + a)-(TA{a)] . (4) 
a 

Here gM corresponds to the hard scattering vertex, M is the heavy particle mass, Xi 
and X2 are the longitudinal momentum fractions of the projectile and target carried 
by the colliding partons, x = xi — X2, = X1X2S where s is the overall energetic 
variable: s = {pi + and p2 is the momentum of a nucleon in the target nucleus. 
Fp(N){.Xi^ M^) is the structure function of the projectile proton (nucleon) at Q"^ = M^. a 
is the light particle- nucleon amplitude, with a = 2 Im a. The notation cr^(a) means the 
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total cross-section calculated with the amplitude a: 

aA{a) = 2 d% (l - [1 + iaTA{b)f) , (5) 

with TA{b) the nuclear preofile function normalized to 1. Finally, we have introduced the 
amplitude a for the scattering of the on-mass-shell heavy system on a nucleon from the 
target nucleus[|, with the cross-section a = 2 Im a. 

An explicit derivation of Eqs. (3) and (4) is given in Sect. 3. Expressions (3)-(4) have 
a simple physical meaning. Note first that for the production on a single nucleon they 
reduce to 

li^'''\x) = ^F,{x„M')F^{x2,M') . (6) 

This expression corresponds to the standard QCD approach in which the heavy flavoured 
system is produced in a hard collision of two gluons. The nuclear effects are just given by 
the total cross-section cr^(a + a) with an amplitude equal to the sum a + a. It corresponds 
to the scattering on the nucleus of a "projectile" containing light and heavy particles. 
The negative term in Eq. (4) eliminates the contribution with no interaction of the 
light particle, which, according to the definition given in the Introduction, is part of the 
internal contribution. As we shall see below the internal contributions have to be treated 
differently and turn out to be numerically small at mid-rapidities. 

In order to see more clearly the physical meaning of Eqs. (3) and (4), we write them 
at fixed impact parameter, assuming that the amplitudes a and a are purely imaginary. 

Eqs. (3) and (4) reduce in this case to 

^In the case when the heavy particle is either a charmonium or bottonium bound state there are 
two different time scales involved, the production of the heavy system containing the cc or bb pair and 
that of the bound state. Although we make no distinction in what follows between heavy system and 
heavy particle it has to be understood that what propagates through the nucleus and interacts with the 
nucleons is not the bound state but the pre-resonance heavy system which, for convenience, we denote 
5' in the following. 
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iTix, b) = ^F,ix„ M')F^{x,, M\ 6)e-^^^^-(^ (7) 

where 



1 r 



1 _ ^~hATA{b) 



has the meaning of the structure function of the target nucleus at impact parameter 
h (see below). The last factor in Eq. (7) corresponds to the rescattering of the heavy 
system. By comparing with the corresponding rescattering in the probabilistic approach, 
Eq. (1), we obtain the substitution (2) mentioned in the Introduction if we neglect the 
shadowing corrections to the nuclear structure functions coming from (8). For J/ip 
production, a ~ aabs and the first absorptive correction (i.e. the term proportional to 
a) is the same in the probabilistic approach and in the asymptotic formula. Since aahs is 
small, the absorptive corrections to J / ip production due to the rescattering of the heavy 
system will be similar in the two cases (see Sect. 4 for a more detailed discussion). This 
result, which was first anticipated in Ref. |]^, is by no means trivial. Indeed if one would 
naively remove the time ordering in the probabilistic expression (performing the integral 
in z from — oo to +oo) one would obtain an absorptive correction of the form (7) but 
with an exponent two times larger. 

Another important point is that (7) contains the total cross-section a instead of the 
absorptive part (Jabs- As a result, the situation is very different for charmonium and for 
open charm production. In the latter case aahs = and thus the probabilistic expression 
gives an behaviour. This result will also be obtained from our field theoretical treat- 
ment in the low energy limit. However, at high energies, as seen from (7), there will be 
absorptive corrections of the same magnitude as for charmonium. This is an interesting 
prediction of our approach. 



In concluding this Subsection we would like to emphasize that the shadowing cor- 
rections to the nucleus structure function Eq. (8), have been described by an eikonal 
formula or, more precisely, by a factorized multi-Pomeron vertex. We could have used in- 
stead the sum of tree diagrams with the triple Pomeron coupling (eikonalized Schwimmer 
formula) or a pure eikonal formula. The smallness of the shadowing corrections in these 
approaches is due to the weakness of the triple Pomeron coupling or to off-mass-shell 
effects (in the latter case). Our treatment can certainly be generalized to include the 
pure eikonal mechanism and, possibly, other mechanisms as well. However our aim in 
this work is not to describe in detail these shadowing corrections but rather the nuclear 
effects due to the rescattering of the heavy system. So we do not elaborate on possible 
generalizations as far as the EMC effect at low x is concerned. This effect is small at low 
energies and vanishes in the low-energy limit. This result, which will be derived explicitly 
in Sect. 3, is due to the fact that at low energies the momentum fraction X2 associated 
to the target is large - outside the shadowing region. Technically, in the formalism of 
Sect. 3, this vanishing is due to tmin effects. 

2.2 The internal component 

We turn next to the internal component. As defined in the Introduction the internal 
component of the projectile (target) corresponds to the case when the heavy system 
is produced from partons belonging to the projectile (target) and no interaction takes 
place between light partons of projectile and target. It turns out that in the central 
region this component is very small at high energies. At low energies it is also small for 
charmonium. For open charm it is not small but is absorbed exactly in the same way 
as the external component. Therefore, in the central region one would obtain the same 
results for the A-dependence by neglecting this internal component. In the projectile 
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fragmentation region, the internal component from the target becomes dominant, since 
other components evidently vanish in the limit xp ^ 1. To calculate it correctly, however, 
one has to properly take into account shadowings corrections to the nuclear structure 
function at extremely low x, which are presumably very important. This will be done in 
our numerical calculations as discussed in Sect. 4. 

The internal contribution evidently involves the probability to find the heavy system 
with a given scaling variable x[ in the projectile - described by the part of the projectile 
structure function contributed by the heavy flavoured particle F^/p(a;'^, Q^)- Since the 
heavy system formed in the projectile may have an arbitrary scaling variable x'l > xi, 
integration over x[ is expected. Upon colliding with the nucleus the heavy system with 
the scaling variable x'l has to transform into the observed heavy system with the scaling 
variable xi. This process is described by the inclusive cross-section 7^^_»$(xi — > xi). 
Prom the diagrams in Fig. 2 one finds 



The heavy particle component of the structure function can be calculated via the light 
particle one assuming that the heavy particle is produced inside the projectile by a hard 
scattering mechanism (i.e. that there is no intrinsic heavy particle component). In our 
simplified scalar case we find 



From (10) one concludes that the heavy flavour structure function is smaller than the 
ordinary one by a factor g"^ (actually taken at the scale ~ and therefore small). 
It is also damped in the region x ~ 1 by an extra power oi {1 — x). 

As to the inclusive cross-section /^^^^(xi — > x), we can estimate it knowing that in 
the rescattering the heavy system tends to conserve its longitudinal momentum. Then 






(10) 
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we approximately have 

/$A->-*(2^1 x) — Xia-q/ASix — Xi) . (11) 

This simple result is however valid only for open heavy flavour {DD, etc.). For hidden 
heavy flavour (e.g. J/ip) one expects an extra absorption in each inelastic collision due 
to transition into open heavy flavour channels. To describe this absorption we introduce 
a factor e which represents the part of the initial heavy flavoured cross-section which 
contains the observed heavy system. For open flavour e — 1 and for charmonium or 
bottonium £ <^ 1. After introducing s, (11) changes to 

A 

hA-*^{xi x) ^ xi5{x - xi) ^ e^a'^l, (12) 

n=0 

where denotes the "if -A cross-section with n inelastic interactions. 

Putting (12) into (9) we obtain the internal contribution from the projectile in a 
simple form (an explicit derivation is given in Sect. 3): 

/r'")(^)=4]ii^*/p(^i,M^), (13) 

where 

= (^'iAi^) + («) - cr*A (« «(1 - ^)) (14) 
and cr^Ai^) is the "^A inelastic cross-section calculated with the "^f-N amphtude d. For 
open flavour production (s = 1) we evidently have = (J^a\ whereas for hidden 
flavour with large absorption {e ~ 0) ~ cr^A- 1^ the last case the internal contribu- 
tion from the projectile is therefore substantially reduced. 

Apart from the internal flavour in the projectile, we have also to take into account 
the internal flavour in the target nucleus, which upon scattering off the projectile will 
also contribute to the cross-section. Its expression can be easily obtained. We get 

4-*'^)(a;) = eaF^/j,{x2, M^)^a, (15) 
10 



where a is the ^ A^" cross-section and the factor is given by (4) . 
For the production on a single nucleon we find 

lf'^^[x)=e~aF^I,{x^,M^) (16) 

and 

7{;"*'^)(a;) = £aF^/^(a;2,M2) . (17) 

To have a more exphcit form we rewrite (16) and (17) for fixed impact parameter h. 
Assuming that a and a are purely imaginary and restricting ourselves to the case of open 
fiavour (e = 1) we have 

lf'^^\x, b) = 2F^/,{x,, M2)(1 - e-(V2)^AT^{b)) (ig) 

and 

I^r''^\x, h) = aF^/A{xu M\ 6)e-(V2)^AT^W^ (19) 

where F^/p(xi, M^, h) is the heavy fiavoured part of the nucleon structure function, Eq. 
(10), at fixed 6, and 

F^/a{x2, M\ h) = F^/^{X2, M2)1[1 - e-(^/^)-^^-(^)] . (20) 

As it was the case for the external component, (18) and (19) actually describe two 
different physical effects: a change of the gluon distribution in the colliding nucleus as 
compared to the free nucleon (shadowing part of the EMC effect) and the rescattering 
of the heavy system in the nucleus target. The former effect corresponds to terms of 
the second and higher powers in the light-particle interaction with the target, i.e. in the 
amplitude a or cross-section a. Rescattering of the heavy system is described by the last 
factor in (18) and (19). 

Note that the internal part coming from the projectile, Eq. (18), is absorbed in a 
different manner. It is proportional to the total cross-section for ^-A scattering and 
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behaves like A^l"^ whereas the two other parts behave hke A^l"^ (at very large A). How- 
ever, as we shall see numerically, both internal contributions are much smaller than the 
external one at Xi? ~ 0. 

In concluding this Subsection we would like to emphasize that, while the internal 
contribution does contain e, the external one does not. This is due to the following. In 
order to compute the inclusive cros-section for the production of the heavy system we 
have to fix an intermediate on-mass-shell heavy system in Fig. 2 (i.e. "cut" the diagram 
through the heavy system line). This cut may pass either through the rescattering blob 
A or through one of the two heavy particle lines with which the blob is attached to 
the rest of the diagram. If the cut passes through the blob A then the lower blob Bt 
may be both cut or uncut. The corresponding contributions cancel |§] in this sum due 
to the well-known AGK cancellations 0. Therefore we have to consider only the two 
cases when the cut passes either through the left or through the right heavy particle 
line in Fig. 2b. Due to this important result the blob A, containing the inclusive cross- 
section of the heavy system, is never cut - and thus e never appears. However, this AGK 
cancellation is only true for the external component. Clearly, it does not take place in 
the case when the light particle does not interact with the nucleus at all. Therefore the 
internal contribution does involve the cut blob A and has to be treated separately. In 
particular it contains the factor e. 

2.3 Generalization to nucleus-nucleus collisions 

In spite of a considerable complication in the dynamics, it can be shown that the AGK 
cancellation which governs the contributions to the inclusive cross-sections discussed in 
the last paragraph of Sect. 2.2 remains valid for AB collisions. Namely, emission of the 
heavy system from the rescattering blob is cancelled out, unless there is no interaction of 
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the light particle with one of the nuclei. As a result, in the AB case we have to consider 
diagrams of the same type as for hA scattering and, correspondingly, the contributions 
to the heavy system inclusive cross-section are divided into an external part, with inter- 
actions of at least one light particle of each nucleus with nucleons of the other nucleus, 
and two internal parts: that of the nucleus B, with no light particle interaction with 
nucleus A, and that of the nucleus A, with no light particle interaction with nucleus B. 
All three parts are calculated quite similarly to the hA case. We present here only the 
final results for the inclusive cross-sections. 

For the external part we find an expression which is a generalization of Eqs. (3) and 
(4): 

Ia"'\x) = ^Fn{xi,M^)Fm{x2,M^)^a^b . (21) 

As one observes, the inclusive cross-section (21) is factorized in the two nuclei, A and B. 
Recalling (3) and (6) it can be rewritten as a simple relation 

j{ext) j{ext) j{ext) 

^AB _ ^B /22^ 

j{ext) j{e.xt) j{ext) ' ^ ' 

J-N -In -'TV 

where in the hA and hB collisions the nuclei are taken in the same kinematics as in 
the AB collision. Evidently (22) means that the total absorption is just the product of 
absorptive factors coming from both nuclei. As is well-known this relation also holds in 
the probabilistic approach 

For the internal parts we also obtain expressions factorized in the two colliding nuclei. 
The internal heavy flavour present in the nucleus B gives a contribution which generalizes 
Eqs.(13) and (14): 

T(int,B)/- ^ __(e) tt, 

' -//N\ 



iTJ'''\^) = 4\f^/n{xi,M')<!>b ■ (23) 



To this internal contribution a similar one J^*^*'^"* has to be added, which takes into 
account the internal heavy flavour of the target. It is given by (23) with the exchanges 
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A ^ B and Xi <-> X2- 



3 Finite energy corrections 
3.1 The formalism 

As discussed in the Introduction, finite energy effects are very important in heavy ffavour 
production. These effects are present up to rather high energies especially at small xp 
(see below). In this Section we are going to consider these corrections. By introducing 
them we obtain equations valid at all energies. At asymptotically high energies they 
coincide with the ones derived in the previous Section. In the low energy limit, they 
coincide with the result of the probabilistic approach, Eq. (1). As discussed in Refs. 
and the finite energy corrections have a clear origin. Due to the presence of the heavy 
system, some of the contributions to the inclusive cross-section have a non-vanishing 
minimal transverse momentum {tmin 7^ 0) and are suppressed by the nuclear form factor 
in a well-defined way. These modified cutting rules have been computed in Ref. ^ in the 
framework of a specific parton model. However, their physical content is so transparent 
that they presumably have a more general validity. They can be summarized in the 
following way. Let us consider a particular ordering of the longitudinal coordinates zj 
of n interactions with the nucleus of the light particle and the heavy system, n = I + h, 
where l{h) is the number of the interactions of the light particle (heavy system) and 

Zi<Z2<---<Zn . (24) 

In this case, at finite energies, the n-th power of the nucleus profile function T4 which 
appears in the expansion in the number of interactions of the cross-section a a, Eq- (4), 

has to be replaced 0] by one of the following integrals: 

/+00 r+00 r+00 
dzi dz2... dznewii^izi- Zj))YlpAib,Zi), (25) 
-00 Jzi Jzn-l i — l 
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where j = 1, 2, n and 

A = niNM^/s x+ . (26) 

Note that for A = 0, corresponding to asymptotic energies, all integrals T^-'^ are equal 
to and are independent of j. 

For A non zero the value of j to be taken in Eq. (25) depends on the particular 
discontinuity of the scattering amplitude we are considering. Namely, all the disconti- 
nuities containing T^^^ are of the following type. Interactions with the nucleus from 1 to 
j — I have to be located to the left of the considered cutting line ("cut to the right"). 
Interaction j may either be cut or be located to the right of the cutting line ( "cut to the 
left"). All the other interactions, from j + 1 to n may be cut in all possible ways. To 
this contribution one has to add its complex conjugate. 

The physical content of these rules is quite clear. When the first interaction, in the 
ordering (24), is cut the exponential damping factor in (25) is not present, i.e. T^^^ = T^. 
Clearly this is the only case where tmin = @]. In all other cases the exponential 
damping factor is present and depends on the longitudinal distance Zj — zi between the 
first interaction 1 and the first one j which is either cut or cut to the left . 

Turning to the inclusive cross-sections, we have to stress that for A > the AGK 
cancellation discussed at the end of Sect. 2.2 is no more valid. Therefore we have to 
consider all possible cuttings on the same footing. 

In treating the emission from the rescattering blob we apply the same approximation 
as was used in Sect. 2 for the internal contribution from the projectile. Namely we 
assume that the observed heavy system conserves the original longitudinal momentum. 
As in Sect. 2 we introduce the factor e < 1 for each inelastic interaction of the heavy 
system to take into account a possible leakeage of the hidden fiavour into the open one. 

Using these simple rules it is now easy to write the expression for the inclusive cross- 
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section at finite energies (i.e. A 7^ 0). In order to do so one has to remember that 
when cutting an interaction to the right (left) the amplitude ia is replaced by ia {—ia*), 
irrespective of whether the particle is light or heavy. Cutting an interaction of a light 
particle one has to replace ia by a. Cutting that of a heavy particle one has to replace 
ia by ea. 

The expression for the external part of the inclusive cross-section can then be written 
in the same form (3), (4), where now 

0-^(0 + a) - (7A(a) — ?A(a, a) — ^^(a = 0, a), (27) 

with 

A n „ 

dAia, a) = Re ^ Q ^ / £bTi^^ {h)a^^ {h) . (28) 

n=l j=l ■' 

Here for j > \ 

a^\h) = 2 (ia + iay-\-ia - ia - {I - e)a)[-{l - £)cF]"-^' (29) 

and for j — 1 

a^\h)^{a + e~a)[-{l-e)~ar' . (30) 

Eq. (29) is quite obvious: the factor [ia + iay~^ results from the cutting to the right 
of the first j — 1 interactions; the second factor corresponds to the cutting of the j-th 
interaction and the last one to the cutting of interactions from j + 1 to n. The first factor 
is trivial since a right cutting does not change the amplitudes ia and ia. The cutting 
of the interaction plus its cutting to the left that appears in the second factor leads to 
—ia for a light particle and to —ia — (1 — £)a for a heavy one. Indeed, the cutting of 
a light particle interaction in all possible ways (right + left + interaction itself) gives 
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ia — ia* + a — 0. The cutting of a heavy system interaction in all possible ways gives 
{s — l)a. This also explains the last factor of (29). 

Eq. (30) is obtained after putting j — 1. In this case the j-th interaction is the first 
one, and thus it has to be cut either to the left or through its interaction. It can be seen 
that the former (cutting to the left) is already included in the complete conjugate term 
responsible for the factor 2 in (29). Only the cut interaction remains in this factor. No 
complex conjugate term appears in this case. 

For the internal part from the projectile, applying the same rules, we obtain Eq. (13) 
where the cross-section is now replaced by a^pA, with 

A n „ 

ai\ = Re E E y d'bTi^\b)a^^^ , (31) 

n=l j=l 

where for j > 1 

= 2iiay-\-ia - (1 - e)a)[-{l - £)a)]"-^- (32) 

and for j — 1 

a^:^ ^ ea[-{l - e)a)r' . (33) 

The internal part coming from the target contains the same rescattering diagrams as 
the external part. This resulted in the same rescattering factor at asymptotic energies, 
(cfr. Eqs. (3) and (15)). Therefore at finite energies it is given by the same Eq. (15) 
with the substitution (27) in 

Since, by construction, each possible discontinuity has been included in one and only 
one term of (28)-(33), it is clear that the above formulae provide an explicit derivation 
of the asymptotic energy results given in Sect. 2. Let us consider the asymptotic case 
when A = and Tp) = TJ. It can be easily checked that one recovers the results of that 
Section. Take Eq. (28). One can notice that, with T^-^^ = T^, there is a cancellation 
between the term of cr*-^-* proportional to —ia — ia and the term of cr'-^"'"^-* proportional 
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to —(1 — £)a. As a result, in the sum over j, one is left with the term proportional to 
—ia — ia from a^'^\ the one proportional to —(1 — s)a from a^'^^ and the term with j = 1. 
The first term gives (7^(0 + 0) and the sum of the other two terms gives —a^J^\{l — e)a). 
Subtracting ^^(ci = 0,5) (see Eq. (27)) one obtains Eqs. (3) and (4). Likewise one can 
show that for A = one recovers Eqs. (13) and (15) for the internal components. 

Finally we shall study the low energy limit by taking A — > 00. In this case only 
terms with j — 1 survive in Eqs. (28) and (31). We obtain from (30), for the external 
contribution, 

^(tLoo)(^) = %F,{x„ M')F,{x„ M^)$5,°) . (34) 
For the internal contribution of the projectile we obtain from (33) 

= ^^F^/Pi^^, M^)$?) (35) 

and for the contribution of the internal heavy flavour of the target nucleus we get from 
(30) 

'(A,A^(X)) 

In all cases 



4Xa^oo)(^)=^^^*/p(^2,M^)$5) . (36) 



^5) ^ ^|!!((1_£M . (37) 
{l-e)a 

Two important observations can be made upon inspection of Eqs. (34)-(37). First, in 
the cross-sections (34) and (36) only terms linear in a (or a) have survived. This means 
that in the low energy limit the screening corrections to the nuclear structure function 
disappear, as expected. Second, comparison of Eqs. (34)- (36) shows that the absorptive 
corrections for all parts of the inclusive cross-section, external and internal, turn out 
to be the same in the low energy limit. This allows to write for the total inclusive 
cross-section on the nucleus in this limit 

(m)//-| \~\ 

^(tLoo)!-) = = ■ (38) 
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Since aabs — (1 — £)cf we recover exactly the probabilistic expression (1). As far as we 
know this is the first time that this expression has been derived in a field theoretical 
approach. 

3.2 Computational methods 

As we see from the above formulae, in order to compute nuclear effects at finite energies 
we have to deal with the new profile functions T^^^ which depend on the longitudinal 
order of the coUisions and involve the parameter A (Eq. (26)). At first sight it seems an 
impossible task due to multiple longitudinal integrations in (25). However the situation 
improves if we take a Fourier transform with respect to A. Then instead of T^^^ we find 
integrals 

Fi^\b,0 = I ^T(^)(&, A)exp(zAO = 

/+0O f+oo r+oo " 

dzi dz2... dzn5{zi - Zj + ^)Y[pA{b,Zi) . (39) 
-oo Jzi Jz„-i j_2^ 

Evidently F^^O = for ^ < and 

Fi'HO-mT'' (40) 

(here and in the following we suppress the argument b and the subindex ^4). For j > 1 
and ^ > simple calculations lead to 

, . 77 1 r+oo 

F!f\b) = dzp{z)p{z + i)T^-\z, z + i)T-'^{z + 0, (41) 

where 

T{z)^ r^dz'piz') (42) 

J z 

and 

T(^i, Z2) = r dz'piz') = Tiz,) - T{z2) . (43) 
Thus the calculation of all nontrivial F^^^ reduces to the one-dimensional integral (41). 
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Moreover, using the representation (41), one can perform the sums over n and i in 
the formulae for the inclusive cross-sections (28) and (31) and thus obtain for them 
closed expressions as a function of ^. We restrict ourselves to the pure rescattering 
contributions, that is to terms linear in a. 

For the external contribution we then obtain 

7j-*) = 4^"*)($J)+X), (44) 

where the first term in the brackets, with given by (37), represents the contribution 
from the term with j = 1 (this is the contribution (34), which survives in the low-energy 
limit) and the second term represents all the terms with j > 1. It is given by an integral 

X = A{A - l)Re j d% j d^e-'^^ dzp{b, z)p{b, z + ^) 

{2id + r)a + idT{b, z,z + i){Aid + {A - l)r)a) - r]a{2id + r]^)T{b, z + C)) (1 + w)^"^ 

(45) 

where we have put 

w = iaT{b, z,z + C)- V^T{b, z + i) (46) 

and T) — 1 — e. 

A similar result holds for the internal contribution from the target nucleus: 

7j"*'^)(x) = 7g"*'^)(x)($J) + X). (47) 
For the internal contribution from the projectile we obtain 

^Jin^t,p)^^m ^y^^ (48) 

where again the first term comes from j = 1 and is the one that survives at low energies 
and the rest comes from j > 1 and is given by an integral 

Y = -2A{A - l)Re ^~a^-^^±^ 

ea 
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J d% J d^e"^^ J ^ dzp{b, z)p{b, 2 + 0(1 + w)"^'^ ■ (49) 

The calculations although rather cumbersome are now feasible. We have taken a 
standard Saxon- Woods form for the nuclear density p(6, z) and replaced all (1 + idT)"^ 
by exp(74zaT). 

4 Numerical results and discussion 

In this Section we present our results for pPh collisions at several energies. For onium 
production we take e = 0.001 and for open charm and bottom e = 0.999. In the 
calculations we take the amplitudes a and d purely imaginary. First, we restrict ourselves 
to the nuclear effects due to the rescattering of the heavy system, i.e. we take only linear 
terms in a or a. 

We start with the asymptotic formulae of Sect. 2. Neglecting the EMC effect at low 

X 

^a = ARa-aJ d%TA{h) exp (^-laAT^(6)) . (50) 

We present our results in terms of Ae// = Ipph/ IpN- The three parts of ^e//, the external, 
internal from the projectile and internal from the target, are given by 



^{ext) _ ARa j^{int,p) _ '^j^l/^'^ j^{int,A) _ ARa 

'^^"l + ri + rs' ~ 1 + 1/ri + r2/ri' ~ 1 + l/ra + n/ra ^ 



where ri and r2 are ratios to the external part of the internal ones contributed by the 

projectile and target, respectively, for a single nucleon as a target. These ratios are 

shown in Fig. 3 for open charm [e = 0.999) taking a =7 mb at s = 60 GeV^ in 

accordance with the data |l^|. For the structure functions Fp jy we have taken the 

GRV LO parametrization |[Tl[] and the structure function Fq,/p has been calculated using 

Eq. (10). One observes that ri is very small for all xp > 0, whereas r2 grows with xp 
^This is the center of mass energy of the collision 'i'-N for experiments at -^3=20 GeV and xf=0- 
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and becomes greater than unity for xp > 0.44-0.5. From (51) we then conclude that the 
internal contribution from the projectile proton is very small for all xp and can safely 
be neglected. The internal contribution from the nuclear target, though also very small 
in the central region, becomes dominant at high xp, since the rest vanishes as xj? — > 1. 
The total Agjf is given by 

Ajf ^ Aijp + 4;'/^ = AR^, (52) 

with the absorptive factor Ra irrespective of the relative weight of the two remaining 
contributions. 

As commented in the Introduction, this absorptive factor is different from the proba- 
bilistic one. However for J/ip production this difference is quite small numerically. One 
can see this from Table 1, where we present A^fj and its three components compared to 
the probabilistic value A^'""^ at xp = and 0.5 and various energies. The dependence of 
a on the energy was taken in accordance with the soft Pomeron picture: a ~ s^"®. 

A completely different result follows for open charm production, for which the prob- 
abilistic approach gives no absorption altogether. Our asymptotic formulas, on the con- 
trary, lead to considerable absorption, as presented in Table 2. However these results 
can only be trusted at very high energies, as we shall presently see. 

To see the finite energy effects we calculated the same quantities using our formalism 
of Section 3. The results are presented in Tables 3 and 4 for J/tp and open charm, 
respectively. Comparing with Tables 1 and 2, we observe that finite energy effects are 
much stronger for the open charm than for the hidden one. Indeed, for open charm at 
low energies the finite energy corrections make the absorption quite small (a ~ 0.98 in 
the A-dependence A°'), bringing the results in accordance with the experimental data [Q. 
This effect only dies out at ^/s > 200 GeV, when our model predicts an absorption for 
the open charm of the same order as for the hidden one. For J ftp the finite energy effects 
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turn out to be rather small. They also reduce absorption, but the effect is insignificant 
(~ 3%). 

As to the Xi?- dependence, in all cases absorption grows with xp due to the growth 
with energy of the ^-N cross-section, although this growth is rather mild. However it is 
strengthened if one takes into account the shadowing corrections to the nuclear structure 
function, which are quite large at very small X2 relevant for heavy flavour production at 
large xp and s. This has been done in a simplified manner multiplying our results by 
the absorption factor due to nuclear corrections to structure functions taken from Refs. 
IIT2I and [1^. The results are shown in the last but one columns of Tables 1 and 2 and 
in the last columns of Tables 3-6. We see that the xp dependence of the EMC effect 
is quite important. Once it is taken into account, the resulting Xi^-dependence of J/ip 
suppression turns out to be consistent with the experimental one [1^. Note, however, 
that at the two lower energies of the Tables, the experimentally observed A-dependence 
of Drell-Yan pair production is very close to at all values of x^. Thus, a detailed 
check of this A-dependence, with our formalism for the EMC effect, is needed before we 
can claim to have a consistent explanation of the Xi?-dependence of J /if) suppression. 
This interesting point is, however, beyond the scope of the present work. 

In Tables 5 and 6 we present the same results for hidden and open bottom production 
respectively taking a = 2 mb at ^/s = 20 GeV and assuming the same energy dependence 
as in the case of charm production. The A-dependence of T at LHC energies and xp=0 
corresponds to a = 0.96, which value is lowered to a = 0.91 when shadowing in the 
nuclear structure functions is introduced. 

Finally results for J ftp production aX xp = in Pb-Pb collisions are presented in 
Table 7. The difference between the asymptotic formula and the probabilistic one is 
larger than for pA but still small. 
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5 Conclusions 



The probabilistic formula used up to now to describe heavy flavour production off nuclei 
has been generalized to all energies using a quantum field theoretical approach. For J/ ip 
and T production it gives practically the same results as the probabilistic formula up 
to i/s ~ 6 TeV. For open heavy flavour production we predict nuclear absorption for 
^/s > 40 GeV. For ^/s >200 GeV this absorption turns out to be almost the same as the 
suppression of the J/ip. Our formalism also predicts an increase of the J/ip suppression 
with increasing xp- 
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Figure captions 



Fig. la. Low-energy heavy flavour production amplitude with rescattering. 



Fig. lb. High-energy heavy flavour production amplitude with rescattering. 



Fig. 2a. Factorized diagram for heavy flavour production without rescattering. 



Fig. 2b. Same as Fig. 2a but with rescattering of the heavy system with partons 
from the nucleus A. 



Fig. 3. Ratios of the internal parts over the external one for projectil (dashed line) 
and target (solid line) contributions for open charm. 
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Tables 



Vs, GeV 


A ext 
^eff 


^eff 


A mt,A 
^eff 


A tot 

^eff 


4tot+EMC 
^eff 


j^prob 


xf—0 


20 


128.9 


3.8 


0.016 


132.7 


144.6 


134.4 


39 


125.7 


2.5 


0.010 


128.2 


134.5 


131.7 


200 


117.6 


1.0 


0.004 


118.6 


101.8 


124.7 


6000 


99.4 


0.3 


0.001 


99.8 


71.0 


109.5 


xf=0.5 


20 


122.3 


0.50 


0.4 


123.2 


117.2 


129.1 


39 


115.9 


0.27 


0.3 


116.5 


96.9 


123.6 


200 


98.6 


0.08 


0.2 


98.9 


70.3 


109.1 


6000 


61.4 


0.01 


0.2 


61.6 


43.1 


78.6 



Table 1: Effective atomic numbers for J/^ production in pPb collisions with asymptotic 
formulae. 
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Vs, GeV 


Aext 


Aint,p 


AintjA 


A tot 


Atot+EMC 


j^prob 




20 


103.1 


16.5 


12.9 


132.5 


144.4 


207.9 


39 


108.3 


11.3 


8.7 


128.3 


131.9 


207.9 


200 


110.6 


4.7 


3.5 


118.8 


101.9 


207.9 


6000 


97.9 


1.2 


0.8 


99.9 


71.0 


207.8 


xp=Q.'o 


20 


28.3 


0.58 


93.9 


122.8 


116.8 


207.9 


39 


32.5 


0.36 


83.4 


116.3 


96.8 


207.9 


200 


28.5 


0.09 


70.4 


98.9 


70.3 


207.8 


6000 


17.1 


0.01 


11.1 


61.6 


13.1 


207.7 



Table 2: Effective atomic numbers for open charm production in pPb collisions with 
asymptotic formulae. 



Vs, GeV 


A ext 
^eff 


tint,p 
^eff 


Aint,A 
^eff 


Atot 

^eff 


Atot+EMC 
^eff 






Xp 


=0 






20 


134.9 


0.1 


0.017 


135.0 


147.2 


39 


133.2 


0.6 


0.011 


133.8 


137.5 


200 


118.4 


1.0 


0.004 


119.4 


102.4 


6000 


99.6 


0.3 


0.001 


99.9 


71.0 






Xp- 


=0.5 






20 


126.4 


0.33 


0.42 


127.2 


121.0 


39 


116.5 


0.26 


0.30 


117.0 


97.3 


200 


98.7 


0.08 


0.24 


99.1 


70.5 


6000 


61.4 


0.01 


0.16 


61.6 


43.1 



Table 3: Finite Energy effective atomic numbers for J/^ production in pPh collisions. 
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Vs, GeV 


Aext 
^eff 


/iint,p 
^eff 


tmt,A 
^eff 


A tot 
^eff 


Atot+EMC 
^eff 






Xp 


=0 






20 


162.7 


20.6 


20.4 


203.7 


222.0 


39 


155.8 


13.5 


12.5 


181.8 


186.9 


200 


114.1 


4.7 


3.6 


122.5 


105.1 


(3000 


97.8 


1.2 


0.8 


99.8 


71.0 








=0.5 






20 


33.8 


0.63 


112.2 


146.7 


139.5 


39 


33.2 


0.36 


85.1 


118.6 


98.7 


200 


28.5 


0.09 


70.3 


98.9 


70.3 


6000 


17.0 


0.01 


AAA 


61.1 


42.7 



Table 4: Finite Energy effective atomic numbers for open charm production in pPh 
collisions. 



^/s, GeV 


A ext 




Aint,A 


A tot 


A tot+EMC 
^eff 






■ft 


=0 






20 


181.5 


0.02 


0.006 


181.5 


197.1 


39 


180.3 


0.06 


0.004 


180.4 


186.5 


200 


175.9 


0.12 


0.002 


176.0 


157.7 


6000 


166.9 


0.04 


0.000 


166.9 


127.5 






xp-- 


=0.5 






20 


178.4 


0.040 


0.2 


178.7 


173.8 


39 


175.1 


0.034 


0.1 


175.3 


153.6 


200 


166.5 


0.012 


0.1 


166.6 


126.7 


6000 


142.5 


0.003 


0.1 


142.6 


82.8 



Table 5: Effective atomic numbers for T production in pPh collisions. 
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Vs, GeV 


Aext 
^eff 


Aint,p 
^eff 


Atnt,A 
^eff 


A tot 

^eff 


Atot+EMC 
^e.f.f 






Xf 


=0 






20 


193.3 


7.0 


6.9 


207.2 


224.4 


39 


192.5 


4.5 


4.4 


201.4 


208.2 


200 


174.4 


1.7 


1.5 


177.6 


159.4 


6000 


166.2 


0.4 


0.3 


167.1 


127.6 






Xp- 


=0.5 






20 


96.5 


0.45 


91.4 


188.3 


186.2 


39 


101.6 


0.26 


74.4 


176.3 


156.6 


200 


97.7 


0.07 


69.0 


166.7 


127.6 


6000 


81.9 


0.01 


60.7 


142.7 


83.4 



Table 6: Effective atomic numbers for open bottom production in pPb collisions. 



Asymptotic Energies 


Vi, GeV 


A ext 
^eff 


A mt,p 
^eff 


A mt,A 
^eff 


Atot 
^eff 




20 


16625.5 


492.8 


2.1 


17120.4 


18078.1 


39 


15809.0 


316.6 


1.3 


16126.9 


17341.2 


200 


13825.1 


123.0 


0.5 


13948.6 


15552.8 


6000 


9896.7 


29.8 


0.1 


9926.6 


11997.9 



Finite Energies 



v/i, GeV 


A ext 
^eff 


A i-nt^p 
^eff 


A •int,A 
^eff 


Atot 
^eff 


20 


18316.1 


13.2 


2.3 


18231.6 


39 


17754.0 


74.2 


1.5 


17829.7 


200 


14014.1 


117.3 


0.5 


14131.9 


6000 


9913.2 


29.9 


0.1 


9943.2 



Table 7: Effective atomic numbers for J/^ production in PhPh collisions at xp—Q- 
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Fig. 3. 
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